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We analyze in detail the problem of gauge invariance of the deeply vitual Compton scattering 
(DVCS) amplitude. Using twist 3 one-gluon exchange diagram contributions and the QCD equa- 
tions of motion, we derive the general gauge invariant expression of the DVCS amplitude on a 
(pseudo)scalar particle (pion, He 4 ). Similarly to the case of deep inelastic scattering, the amplitude 
does not depend on the twist-3 quark-gluon correlations at the Born level. The contribution of the 
derived amplitude to the single-spin asymmetry with longitudinally polarized lepton is calculated. 



Deeply Virtual Compton Scattering (DVCS) has re- 
cently attracted much attention. One of the main rea- 
sons of this interest is the fact that the DVCS process 
gives information about a new type of parton distribu- 
tions, called skewed parton distribution (see for example 
El-B and references therein). The process 



(1) 



has been shown to factorize in the Bjorken region with 
(q 1 ) 2 — 0, — q 2 large and small transfer t = (p — p') 2 , 
as the product of a perturbatively calculable coefficient 
function and a long distance object, the skewed parton 
distribution, which generalizes the notion of parton dis- 
tributions. 

The fact that there is a problem with the photon gauge 
invariance of the DVCS amplitude in leading order at 
Bjorken limit is fairly well-known (see, for instance 0). 
The relevant terms are proportional to the transverse 
component of the momentum transfer and provide the 
leading contribution to some observables, and in partic- 
ular, to the Single Spin Asymmetry. 

As was shown in ^j, a fruitful analogy between the 
transverse spin case of the deep inelastic scattering (DIS) 
and the DVCS process can be used to derive a general 
solution of this problem. We elaborate on this approach 
in the current paper. 

For simplicity, we concentrate here on the DVCS pro- 
cess off (pseudo)scalar hadrons, which may be pions or 
hclium-4 nuclei, but our calculation may be generalized 
to any hadrons. We also neglect hadron masses effects, 
i.e. kinematical power corrections, which may be studied 
independently. 

The Lorentz structure of the hard subgraph of the lead- 
ing order DVCS diagram (Fig. la) has the form of a trans- 
verse projector g» u - P^n v - n^P" J|. All 4- vectors 
may be presented in the form of the Sudakov decompo- 
sition over two light-cone vectors P, n and one compo- 
nent transversal to the given light-cone vectors. Hence, 
if the virtual photon momentum, which has a transverse 



component, is convoluted with the hard part of the lead- 
ing order DVCS amplitude, one obtains a term directly 
proportional to the transverse component of the virtual 
photon momentum. In other words, the measure of the 
photon gauge invariance violation is the non-zero trans- 
verse component of the virtual photon momentum. The 
recent general analysis || confirmed that violating terms 
arc indeed kincmatically subleading. 

In this paper we generalize the Ellis-Furmanski- 
Petronzio (EFP) factorization scheme to the non- 
forward case and calculate the complete expression for 
the DVCS amplitude up to twist 3 order. While this 
scheme was originally applied to the case of the twist- 
4 power corrections to DIS, our analysis will be more 
close to the subsequent treatment of the transverse po- 
larization in DIS at twist 3 level ||||, Analogously to 
that case, the process (Q) is described by the diagrams 
of Fig.l, which contain leading and next to leading twist 
contributions. Namely, the diagram (b) , in the case of 
the transverse gluon field, is entirely at twist 3 level, while 
the handbag diagram (a) contains, besides the standard 
twist- 2 term (produced by the good component of quark 
fields and collinear parton momenta), a twist-3 term, re- 
lated to the quark gluon contribution of diagram (b) by 
the equations of motion. The latter play a crucial role in 
DIS, guaranteeing that the sum of diagrams (a) and (b) 
is gauge invariant and, moreover, not explicitely depen- 
dent of the quark-gluon correlations. We will show that 
the situation in DVCS is quite similar. 

The sum of T$ amplitude from diagram (a) and T$) 
from diagram (b), has the following form, like in 
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T (b) = 



dk 1 dk 2 tr\ E f _ ipu (k ll k2)T p (ki,k 2 ) 



(2) 



where E^ u and E^py are the coefficient functions with two 
quark legs and two quark and one gluon legs, respectively. 
For simplicity, we restrict to the Born diagrams for the 
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coefficient functions. In (g), the following notations are 
introduced 



T a p{k) =-j dze t( - k --*> z ( P '\i; a (z)M0)\p), 

r^(fci,fc 2 ) = - J dz 1 dz 2 e i ^-^ +i ^- k ^ 

{p'\Mzi)gA p (z 2 )M0)\p), 
A A 

p' = P+j, p = P-~, A = q-q' 



(3) 



Here, p' and p are the final and initial hadron momenta, 
q' and q are the final and initial photon momenta. For the 
sake of convenience, we neglect all the kinematical power 
corrections and put P 2 and t = A 2 both equal to zero, 
keeping only the terms linear in Ax. Also, we choose the 
axial gauge condition for gluons, i.e. n-A = 0, where n is 
a light-cone vector, normalized by the condition n-P = 1. 
It is convenient to assume that n = q' / P ■ q' , although 
our result will not depend explicitely on the choice of n. 
We carry out a decomposition of k in the basis defined 
by the P- and n- light-cone vectors 

k = X P + (k ■ P)n + k T) X = k-n. (4) 

Apart from this, 

A = -2£P + A T . - 2£ = A • n. (5) 

Further, we carry out the following replacement for the 
integration momentum in eqn. (0) 



dki 



dkidxt8(xi — hi ■ n). 



(6) 



Expanding the two-quark coefficient function E pu (see, 
(||)) in a Taylor series and, next, using the following Ward 
identity 



dkP 



(k,k) 



(7) 



we can write the DVCS amplitude as 
T$+T$ = J dxtr^E^(xP)T(x) 
J dx 1 dx 2 tr^E pp „(x 1 P,x 2 P)uj ppl T p ,(x 1 ,x 2 )^ (8) 
where u! pp > = 5 pp i — n p >P p , and 

T aP (x) = - J dAe^+«Vl^(An)#(0)b), 
T p ^{ Xl ,x 2 ) = - J d\,dX 2 e ^i+mi+<(*2-*i)A 2 
(p'\^ a (X 1 n)D"'(X 2 n)^(0)\p), (9) 

where D p is the QCD covariant derivative in the funda- 
mental representation. 



Let us now focus on the QCD equations of motion 
both for incoming and outcoming quarks. From these 
equations, we deduce integral relations for structure func- 
tions, parametrizing quark and quark-gluon correlations. 
So, let us to start from the QCD equations of motion (we 
consider massless quarks) 



(D(z)1>(z)#(0)) = M*r$(0)£>(0)) = 0, (10) 

where (. . .) denote the asymmetrical matrix elements. 
Keeping only vector and axial projections (since for mass- 
less quarks all other structures do not contribute), we 
decompose the quark and quark-gluon correlators in the 
7-basis. We have 

-4<V>(^(0)) = (4>(0h a iP(z)) la - (^(0)7 Q75 ^(z))7a75 

(11) 

-4{gA"{y)1>(z)$(p)) = (^(0) la gAf{y)^(z)) la - 

W0)7a755^ P (y)V^)>7a75- (12) 

In terms of P,A T ,n- vectors, we introduce the 
parametrization of relevant vector and axial correlators 
(see, Eq. (^Tj, [TJI)) in the following forms, where terms 
proportional to n, contributing only at twist-4 level, have 
been omitted, and the axial gauge condition n- A = has 
been taken into account : 



(tK0)7^(z)} = ffi(z)P M + H 3 (x)Ai 



(m^dj^z)) £ iHj{x)P p A T p 



(7/>(0)757 M VM) = iH A {.x)e llA Tp n 



(^(0)757/^' H z )) = -H'A{x)P p e pA T Pn , 



(mi^iyMz)) £ B{ Xl ,x 2 )P^ p , 



(ip(0)j 5 j p gA p (y)ip(z)) = iD{xi,x 2 )P p e pA T 



Fm 



(13) 
(14) 
(15) 
(16) 
(17) 
(18) 



here £ p A T p n = Spap-y A Ta P l3 n 1 ;= denotes the Fourier 
transformation with measure (z = An, z' = 0) 

dxe' i{xP ~^ )z+l(xP+ ^ )z ' 
for quark correlators, and 

dxidx 2 e~ i< - XlP ~^ s>z ~ i< - X2 ~ Xl ^ Pv+i ( X2P+ ^" ,z ' 

for quark-gluon correlators. Note that the latter for the 
non-forward case are actually new objects. We kept 
only the argument x for all the correlators, dropping 
for brevity the dependence of the distributions on the 
skewedness parameter £, recovering it below where it is 
necessary. Their dependence on t — A 2 is beyond our 
scope. Finally, their dependence on the factorization 
scale parameter u 2 requires an extended separate inves- 
tigation. 

The dependence on £ plays a crucial role in the disap- 
pearance of the i?3 term, when a local current, related to 
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its integral in x, is considered. Although this is required 
by the conservation of local vector current, this effect is 
proportional to t and is therefore beyond the scope of our 
approximation. At the same time, this is also required by 
T-invariance [13]. ft is therefore natural, that symmetry 
of H (c.f. [0,11), resulting also from T-invariance, 

H 3 (x,0 = -H 3 (x,-0 

is relevant. Calculating the integral in x and implying the 
polynomiality condition ||, one gets a function which is 
independent on £, and hence vanishes since the only odd 
constant is zero. A similar argument is also applicable to 
the function H A , so that : 



J dxH 3 (x) = , J dxH A {x) 



= 0. 



Acting id on (|ll|) and 7 on (|12| ) from left or right side, 
the QCD equations of motion yield the following integral 
relations for structure functions 

dy (B^ (x, y) - D^( Xl y) + 8{x - y)H T A {y) 
-ZHzW-^H^-xHAix), 

dy (b^ (a:, y) - 5(x - y)H?(y) - (x, y)) = 

xH 3 (x)+£H A (x), (20) 

where symmetrical and anti-symmetrical functions are 
defined as, 

B^ A \x,y) = 1 (B(x,y)±B(y,x)) . (21) 

Note again the important difference with DIS, where 
the axial correlator is symmetric and the vector one is 
antisymmetric The latter property is based on T- 
invariance, just like the symmetry properties in £ dis- 
cussed above. To see the relation between x <-> y and 
£ symmetry, it is instructive to write the general T- 
invariance relations : 

B(x, y, = B(y, x, -£), D(x, y, £) = -D(y, x, -£); (22) 

So the "unnatural" symmetry in x, y results from the an- 
tisymmetrical in £ part, clearly absent in the forward 
case. It is worthy to note that a similar unnatural 
symmetry may appear due to the final state interaction 
phases in the case of T-odd fragmentation functions p2| . 

We can thus write the DVCS amplitude in a gauge 
invariant manner. Let us first write the contribution from 
the pure quark amplitude : 

1 



rp(a) _ 
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dx- 



where the following notations were introduced 

Sfl 1 fl 2 fJ. 3 fJ4 ~ 9fJ\fJl9fl3fJ.4 + 9 fJ.lfJti9 1*21*3 ~9fJ.lfJ.s9fi2lJ.41 

Q=(q + q')/2 

As a corollary, the contribution of the amplitude cor- 
responding to the one-gluon exchange diagram, has the 
form 



T^ = \Jdx ld x 2 — 



1 



v P + Q) 2 {x 2 P + Q) 2 

V>,»)-s(*>-^M 



(19) ti[-y 1J {x 2 P + Q)A T (x 1 P + Q)-y ll P 



i[ D{xi,x 2 ) - 5(xi - x 2 )H A {x 2 ) )e aA T 



tv[^(x 2 P + Q) la { Xl P + Q^^PjrA )+"crossed". (24) 



Calculating all traces in (|24|), using the following obvious 
identities 

±(P ■ Q){x 1 + x 2 ) + Q 2 



(x 1 P±Q) 2 {x 2 P±Q) 2 
1 { 1 1 



2\ v (xiP±Q) 2 (x 2 P±Q) 2 

±{P ■ Q){ Xl - x 2 ) = 
( Xl P±Q) 2 (x 2 P±Q) 2 

l(l 1 



2\{x 2 P±Q) 2 { Xl P±Q) 2 )' 

and the equations of motion (|2^) in terms of symmetric 
and anti-symmetric functions, we add contributions of 
(^||) and (|24|), taking into account the crossed diagrams. 
The gauge invariant expression of the DVCS amplitude 
is thus 



where 



1 



2P-Q 



dx 



1 



1 



x — £ + ie x + £ — it 



(25) 



xP + Q) 2 



= H x (x) -2£P M P„ - P^Q, - P V Q^ 



g^{P-Q)--PAl^- 2 P^l 



Hl{x)Sv(xP+Q)fjP + I jT 3,{x)S u ( x p + Q) llA T + 



H A (x)e aA r Pn e u ^ xP+ Q) tia + (fi — > v, Q — > -Q) , (23) 



^H A (x) 1 3^Al - £P u Al - AlQ v + A T V Q, 
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We can see that the first term of (p5|), proportional to 
the iii-function, completely coincides with the improved 
DVCS amplitude, proposed by Guichon and Vanderhae- 
gen (GV) in M. Indeed, if we proceed, in such terms, to 
GV basis Q , where the n- vector is expressed via the vir- 
tual photon momentum q and decompose the Q- vector in 
this basis, then we derive, omitting the terms with H3- 
and iHa- functions : 



6( x ' + 0-6(x'-0 



V V 



x — £ x + £ 
H x {x)H 3 {x') - H 1 {x')H 3 {x)+ 



(28) 



T v = S dx \ 



1 



( 

Hi{x) I P^riy + P v n^ - g^ v - — ^ 



T, 



L.O. 



P, 



P-q 



V \TrpL.O. 



(26) 



where the definition of T^ v °- 



amplitude is related with 
the transverse direction projector, which is {P^n v + 

In complete analogy to DIS, the answer does not de- 
pend explicitely on quark-gluon correlations. However, 
in contrast to DIS, it contains two additional new func- 
tions, instead of the single function gi in the DIS case. 
While H3 may be considered as an analog of gi (the co- 
efficients of them are At and st, respectively), the DIS 
analog of the function Ha is excluded by T-invariancc 
and may be present only for fragmentation. 

We emphasize that the deduced gauge invariant ex- 
pression of the DVCS amplitude has a significant mean- 
ing for the sequential application of QCD to the investi- 
gation of any observable values. To demonstrate this, let 
us consider the single (electron) spin asymmetry (SSA), 
which arises in the collision of the longitudinaly polarized 
electron beams with the unpolarized scalar target. The 
SSA parameter (for details see, for instance, 0) is 



where 
da(~ 



At 



da(-+) - da(<-) 



da(- 



dcr(< 



e 6 F+(t)2£ 



q 2 t{k~ A) 2 (fc' + A) : 



;£kk'PA 



dx[S(x + Z)-6(x-0\ ■ [ Hx(x)((k + tf) • P)+ 



2H 3 (x) U((k + k') ■ P) + (k' ■ A) - £(P • Q)\ + 



2i 



x(P-Q) 



H A (x) (k ■ A)(k' ■ P) - (k' ■ A)(k • P) + 



2i 



4 £fefc ' PA ( P . Q) 



dxdx 1 



S(x + 0- S(x - 



V 



V 



x' -i x 1 + £ 



where F+(t) is the target electromagnetic form factor, 
emanating from the Bethe-Heitler diagrams, k and k' 
denote the momenta of the initial and final electron 
(k-k' = q). 
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Fig.l The DVCS diagrams ( notations: k = xP - A/2, 
k' = xP + A/2, m = xiP - A/2 and I = x 2 P + A/2). 
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